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Abstract 



Using the CLEO-II data set we have searched for the Cabibbo-suppressed 
decays /)(*)-. For the decay D*+ D*~ , we observe 

one candidate signal event, with an expected background of 0.022 ib 0.011 
events. This yield corresponds to a branching fraction of B{B^ ^ D*^ D*~) = 
(5.3^1:3, (stat) ± l.O(syst)) x lO""^ and an upper limit of B{B^^D*+ D* ) < 
2.2 X 10-3 at the 90% CL. For S° D*^ and B^ , no sig- 

nificant excess of signal above the expected background level is seen, and 
we calculate the 90% CL upper limits on the branching fractions to be 
B{B^^D*^D^) < 1.8 X 10-3 and B{B°-^D+ D') < 1.2 x 10-^. 



1 



Typeset using REVT^ 



D. M. Asner/ D. W. Bliss/ W. S. Brower/ G. Masek/ H. P. Paar/ V. Sharma/ 
J. Gronberg,^ R. Kutschke,^ D. J. Lange,^ S. Menary,^ R. J. Morrison,^ H. N. Nelson,^ 
T. K. Nelson,2 C. Qiao,^ J. D. Richman,^ D. Roberts,^ A. Ryd,^ M. S. Witherell,^ 
R. Balest,3 B. H. Behrens,^ K. Cho,^ W. T. Ford,^ H. Park,^ P. Rankin,^ J. Roy,^ 
J. G. Smith, ^ J. P. Alexander,'^ C. Bebek,^ B. E. Berger,^ K. Berkelman,'^ K. Bloom,^ 
D. G. Cassel,^ H. A. Cho,^ D. M. Coffman,^ D. S. Crowcroft,^ M. Dickson,^ P. S. Drell,^ 
K. M. Ecklund,^ R. Ehrlich,^ R. Elia,^ A. D. Poland,^ P. Gaidarev,^ B. Gittelman,^ 
S. W. Gray,^ D. L. Hartill,^ B. K. Heltsley,^ P. I. Hopman,^ J. Kandaswamy,^ 
N. Katayama,^ P. C. Kim,"^ D. L. Kreinick,^ T. Lee,^ Y. Liu,"^ G. S. Ludwig,^ J. Masui,"^ 
J. Mevissen,^ N. B. Mistry,^ C. R. Ng,^ E. Nordberg,^ M. Ogg,^'[| J. R. Patterson,^ 
D. Peterson,^ D. Riley,^ A. Soffer,^ C. Ward,^ M. Athanas,^ P. Avery,^ C. D. Jones,^ 
M. Lohner,^ C. Prescott,^ S. Yang,^ J. Yelton,^ J. Zheng,^ G. Brandenburg,^ R. A. Briere,' 

Y.S. Gao,^ D. Y.-J. Kim,^ R. Wilson,^ H. Yamamoto,^ T. E. Browder,^ F. Li,^ Y. Li,^ 
J. L. Rodriguez,^ T. Bergfeld,^ B. I. Eisenstein,^ J. Ernst,^ G. E. Gladding,^ G. D. Gollin, 
R. M. Hans,^ E. Johnson,^ I. Karliner,^ M. A. Marsh,^ M. Palmer,^ M. Selen,^ 
J. J. Thaler,^ K. W. Edwards,^ A. Bellerive,^° R. Janicek,^^ D. B. MacFarlane,^° 
K. W. McLean,^o p Patel,^° A. J. Sadoff,^^ R. Ammar,^^ p Baringer,^^ a. Bean,!^ 

D. Besson,^2 p) Coppage,^^ C. Darling,^^ Davis,^^ N. Hancock,^^ g Kotov,^^ 
1. Kravchenko,^^ N. Kwak,^^ g_ Anderson,^^ y. Kubota,^^ M. Lattery,^^ j_ j_ O'Neill,^^ 
S. Patton,^3 Poling,i3 T. Riehle,^^ y Savinov,!^ A. Smith,^^ g Alam,i^ 
S. B. Athar,!^ Z. Ling,^^ A. H. Mahmood,^^ H. Severini,^^ S. Timm,!^ F. Wappler,!^ 
A. Anastassov,^^ S. Blinov,^^'! J. E. Duboscq,^^ K. D. Fisher,^^ D. Fujino,^^'| R. Fulton,^^ 
K. K. Gan,i5 T. Hart,^^ K. Honscheid,^^ H. Kagan,!^ R. Kass,!^ J. Lee,^^ M. B. Spencer,^' 

M. Sung,i5 A. Undrus,i5't R. Wanke,!^ A. Wolf,^^ M. M. Zoeller,!^ B. Nemati,^^ 
S. J. Richichi,i6 W. R. Ross,^'^ P. Skubic,^'^ M. Wood,!^ M. Bishai,^^ J. Fast,^^ E. Gerndt,^ 
J. W. Hinson,^^ N. Menon,^^ D. H. Miller,i^ E. I. Shibata,^^ I. P. J. Shipsey,^^ M. Yurko,^' 
L. Gibbons,!^ S. D. Johnson,!^ Y. Kwon,!^ g_ Roberts,^^ E. H. Thorndike,!^ C. P. Jessop,^' 
K. Lingel,^^ H. Marsiske,^^ M. L. Perl,^^ g p gchaffner,!^ D. Ugolini,!^ R. Wang,!^ 
X. Zhou,i9 T. E. Coan,20 V. Fadeyev,^^ I. Korolkov,^^ Y. Maravin,^^ I. Narsky,^^ 
V. Shelkov,20 J. Staeck,^^ R. Stroynowski,^^ I. Volobouev,^^ J. Ye,^^ M. Artuso,^! 
A. Efimov,2i F. Frasconi,2i M. Gao,^! M. Goldberg,^! D. He,^! S. Kopp,^! G. C. Moneti,^! 
R. Mountain,2i Y. Mukhin,^! S. Schuh,^! T. Skwarnicki,^! S. Stone,^! G. Viehhauser,^! 
X. Xing,2i J. Bartelt,22 S. E. Csorna,^^ V. Jain,^^ S. Marka,^^ A. Freyberger,^^ 
R. Godang,^^ K. Kinoshita,^^ I. C. Lai,^^ P. Pomianowski,^^ S. Schrenk,^^ G. Bonvicini,^^ 

D. Cinabro,24 R. Greene,^^ L. P. Perera,^^ B. Barish,^^ M. Chadha,^^ S. Chan,^^ 
G. Eigen,25 J. S. Miller,^^ C. O'Grady,^^ M. Schmidtler,^^ J. Urheim,^^ A. J. Weinstein,^^ 

and F. Wiirthwein^^ 

^University of California, San Diego, La Jolla, California 92093 



*Permanent address: University of Texas, Austin TX 78712 
^Permanent address: BINP, RU-630090 Novosibirsk, Russia. 

•^Permanent address: Lawrence Livermore National Laboratory, Livermore, CA 94551. 



2 



^University of California, Santa Barbara, California 93106 
^University of Colorado, Boulder, Colorado 80309-0390 
^Cornell University, Ithaca, New York 14853 
^University of Florida, Gainesville, Florida 32611 
^Harvard University, Cambridge, Massachusetts 02138 
'''University of Hawaii at Manoa, Honolulu, Hawaii 96822 
^University of Illinois, Champaign-Urbana, Illinois 61801 
^Carleton University, Ottawa, Ontario, Canada KIS 5B6 

and the Institute of Particle Physics, Canada 
lOMcGill University, Montreal, Quebec, Canada H3A 2T8 
and the Institute of Particle Physics, Canada 
^^Ithaca College, Ithaca, New York 14850 
^^University of Kansas, Lawrence, Kansas 66045 
^^University of Minnesota, Minneapolis, Minnesota 55455 
^^State University of New York at Albany, Albany, New York 12222 
^^Ohio State University, Columbus, Ohio 43210 
^^University of Oklahoma, Norman, Oklahoma 73019 
^^Purdue University, West Lafayette, Indiana 47907 
^^University of Rochester, Rochester, New York 14627 
^^Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309 
2° Southern Methodist University, Dallas, Texas 75275 
Syracuse University, Syracuse, New York 13244 
^^Vanderbilt University, Nashville, Tennessee 37235 
Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
^^Wayne State University, Detroit, Michigan 48202 
California Institute of Technology, Pasadena, California 91125 



3 



The decays B^^D^*^+ DW" are favorable modes for studying CP violation in B decays. 
In the Standard Model, time- dependent asymmetries in the decays can be related to the angle 
(3 of the unitarity triangle [|l| . This angle can also be measured with K'g decays; any 

difference between the values obtained in Z)*^*^~ decays and B^ ^ ip Kg would 

indicate non-Standard Model mechanisms for CP violation [QJ^. Although B^ ^ D*^ D*~ 
and B^ ^D*^ are not pure CP eigenstates, estimates indicate that a dilution of the CP 
asymmetry of only a few percent would be incurred by treating these modes as pure CP 
eigenstates [0]. 

The modes B^ ^ D^*^+ D^*^' have never been observed, and no published limits on their 
branching fractions exist. The decay amplitude is dominated by a spectator diagram with 
b — s> cW^ followed by the Cabibbo-suppressed process cd. One can estimate the 

branching fractions for 5° £)(*)+ by relating them to the Cabibbo-favored decays 

}3{B'^^D'^*^+ D^*^-) ^ tan^^c (1) 

where the fx are decay constants and 6c is the Cabibbo angle. Table | shows the 
expected B^ — > branching fractions, where the CLEO measurements of 

have been used §. 



TABLE I. Estimated branching fractions for ^ D^*^^ D^*^ based on the measured branch- 
ing fractions of the Cabibbo-favored decays oi*^^ D^*^^. 



Mode 


B of Related 


Estimated B for 




Mode (%) 


d(*)+d(*^- (10-4) 




2.4 


9.7 




2.0 


8.1 


B^^D+ D- 


1.1 


4.5 



The data used in this analysis were recorded with the CLEO-II detector |^ located at 
the Cornell Electron Storage Ring (CESR). An integrated luminosity of 3.09 fb~^ was taken 
at the T(4S) resonance, corresponding to approximately 3.3 x 10^ BB pairs produced. 

At the T(4S), the BB pairs are produced nearly at rest, resulting in a spherical event 
topology. In contrast, non-BB, continuum events have a more jet-like topology. To select 
spherical events we required that the ratio R2 of the second and zeroth Fox- Wolfram moments 
@] be less than 0.25. 

We required charged tracks to be of good quality and consistent with coming from the 
interaction point in both the r — cj) and r — z planes. We defined photon candidates as 
isolated clusters in the Csl calorimeter with energy greater than 30 MeV in the central 
region (cos 6' < 0.71, where 9 is measured from the beamline) and greater than 50 MeV 
elsewhere. Pairs of photons with measured invariant masses within 2.5 standard deviations 
of the nominal vr" mass were used to form vr*^ candidates. Selected vr^ candidates were then 
kinematically fitted to the nominal 7r° mass. 
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A particle identification system consisting of dE/dx and time-of-fliglit was used to dis- 
tinguisli cliarged kaons from charged pions. For charged pion candidates, we required the 
hkehhood of the pion hypothesis, L^r, to be greater than 0.05. Since all signal modes re- 
quire two charged kaons, the kaon candidates were required to have a joint kaon hypothesis 
likelihood, LkiLk2, greater than 0.10. 

We reconstructed all D*^ candidates in the mode D*~^ ^tt~^ (charge-conjugate modes 
are implied). candidates were reconstructed in the modes D^^K~7t^, D^^K~tt~^7:^ 
and — »• K~ tt"'" tt" tt"'". candidates were reconstructed via —>■ 7r+ 7r+. Table |I| 
summarizes the branching fractions of the D^*^ modes used 

TABLE II. Branching fractions of D^*^ modes used in reconstruction. 



Decay Mode Branching Fraction (%) 









68.3 ± 1.4 


D^- 




7r+ 


3.83 ±0.12 


D^- 






13.9 ±0.9 






7r+ vr^ 7r+ 


7.5 ± 0.4 



D+'^K--K+'K+ 9.1 ±0.6 



For the decay mode K~ 7r+ 7r°, we make a cut on the weight in the Dalitz plot in 

order to take advantage of the resonant substructure present in the decay. The cut choosen 
was 76% efficient for good 7r+ vr^ decays while rejecting 69% of the background. 

We performed a vertex-constrained fit on all the charged tracks in the candidate 
for modes that contained a D*~^. The from the vertex fit was required to be less than 
100. The fit improved the determination of the angular track parameters for the slow tt^ 
from the D*~^ decay. The resulting r.m.s. resolution on the reconstructed mass difference 
Amo*-D = m£,*+ — m£,o was approximately 0.69 MeV. 

Because 5° D*^ is a Pseudoscalar — > Vector + Pseudoscalar decay, the cosine 
of the decay angle, cos^7r+, of the slow tt"*" from the D*^ has a cos^^ distribution, while 
background events have a uniform distribution in this variable. For —^D*^ candidates 
we required | cos 6*^+1 > 0.5. 

To select candidates that contain well-identified D^*^s we combine the reconstructed 
D^*^ masses into a single quantity, Xm- The definition of Xm for each mode is given by 

^UD--D'-) = ((A-").-(A'")y4(Am).-(A,n)y^ 

/ \ 2 / \ 2 



+ r " " " (2) 



^15+ / \ ^^0 + 
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where the values in angle brackets represent the nominal values and the sigmas are the r.m.s. 
resolutions on the given quantity. We require x\i{D*^D*~) < 8.0, Xm{D*^ D"^) < 4.0 and 
x\i{D^ D~) < 2.0. From studies of Monte Carlo and regions in the data outside of the signal 
areas in other variables, we find that the backgrounds are uniform in xti- 

Since the energy of the is equal to the beam energy at CESR, we used the beam energy 
instead of the measured energy of the candidate to calculate the beam-constrained mass: 



rriB = yE'^eam ~ Pb- The r.m.s. resolution in for signal events, as determined from 
Monte Carlo, is 2.8 MeV. In addition, the energy difference, AE = Eb — Ebeam, where Eb 
is the measured B^ energy, was used to distinguish signal from background. The resolution 
in AE is 12 MeV after performing a mass-constrained fit that included the masses of all 
secondary particles {D'^*^ and 7r°). The signal region in all modes was defined as \AE\ < 2a ae 
and \mB - {mBo) \ < 2(7^3- 

We used a Monte Carlo simulation of the CLEO-II detector to optimize all cuts. Since 
the number of observed signal events was expected to be small, all cuts were optimized to 
minimize the probability that the expected background level would fluctuate up to or beyond 
the expected signal level. For calculating the expected number of signal events during this 
optimization we assumed a branching fraction of 0.1% for all 5° ^ modes. 

Using the cuts defined above, we determined the signal reconstruction efficiency using 
Monte Carlo. The reconstruction efficiency and single event sensitivity {SES = {e B Nbb)~^, 
where e is the detection efficiency, B is the product of the daughter branching fractions and 
Nbb is the number of BB pairs produced in the data set) for each mode are summarized 
in Table The systematic uncertainty on the SES is dominated largely by uncertainties 
in the D and D* branching fractions, and, due to the large mean multiplicity of the final 
states, the uncertainty in the tracking efficiency. 



TABLE III. Summary of reconstruction efficiencies and single event sensitivities for the three 
50^2?{*)+^{*)- modes. 



Mode 


Efficiency, e 


SES ={eB NbbY^ 




(%) 


(10-4) 


B^^D*+D*- 


1.86 


5.45 ±0.99 


B0^B*± j^t 


5.07 


3.79 ± 0.53 


B^^D+ D- 


14.41 


2.52 ± 0.40 



The dominant background is due to random combinations from BB and continuum 
events. The Monte Carlo predicts that this background varies smoothly in AE and niB , 
and this is verified in the data. The tub distribution for data in AE sidebands (50 MeV < 
\AE\ < 400 MeV) varies smoothly with no peaking in the signal region. The same is true for 
the AE distribution for data with < 5.27 GeV. To estimate the background in the signal 
region, we count the events in a sideband in the AE-itlb plane (50 MeV < \AE\ < 400 MeV; 
5.2 GeV < niB < Eheam) and multiply by the relative efficiencies of the signal and sideband 
regions determined from background Monte Carlo. 

Figures |l|, ^ and ^ show the resulting plots of AE vs. ttib for the three modes. The signal 
region is indicated with a solid line, and the sideband region is indicated with a dotted line. 
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Table |V| lists the event yields in the sideband and signal regions. The expected number 
of background events in the signal region is also given. The uncertainty on the expected 
number of background events is a combination of statistical error on the number of events 
in the AE sideband regions and the uncertainty in the background shape through the signal 
region. 



TABLE IV. Summary of events found in the data, both in the Ai? sidebands and in the signal 
region, for each of the three modes. 



Mode 


Events in AE' 
Sidebands 


Predicted Background 
in the Signal Region 


Events found in 
Signal Region 




4 


0.022 ±0.011 


1 


B0^B)*± B)T 


117 


0.64 ±0.10 


2 




539 


2.64 ±0.34 


3 



The probability that the expected background of 0.022 ± 0.011 events in B^^D*^ D*~ 
fluctuates up to one or more events is 2.2%. If we interpret the one observed event as 
evidence for a signal, the resulting branching fraction would be 

^ D*+ D*-) = (5.3+^:^(stat) ± l.O(syst)) x 10'\ (5) 

where the systematic uncertainty comes from the uncertainty in the SES. 

No significant excess of events is seen in the other two modes. We calculate upper limits 
on the branching fractions for all three modes, and these results are summarized in Table 0. 
The systematic uncertainty in the SES and the uncertainty in the background level have 
been incorporated into the upper limits [^. 



TABLE V. Summary of upper limits on the - 


^ branching fractions. All upper 


limits are quoted at the 90% confidence level. 




Mode 


Upper Limit (90% CL) 


B0^D*+ D*- 


2.2 X 10-^ 


B0^D*± 


1.8 X 10~3 


B^-^D+ 


1.2 X 10^3 



We have performed a search for the decays D^*^^ . In the mode ^ D*^ D*^, 

one event is seen in the signal region where the expected background is 0.022 ± 0.011. The 
one event in B^ D*^ D*~ is seen at a rate that is consistent with predictions, and in all 
three modes the upper limits are within about a factor of two from the predicted branching 
fractions. 
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FIG. 1. AE vs. niB for data in the ~^D*^ D* analysis. The signal region is indicated by a 
solid box. The sideband region lies above the top and below the bottom dotted lines. 
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FIG. 2. AE vs. niB for data in the ^D*^ D^analysis. The signal region is indicated by a 
soUd box. The sideband region hes above the top and below the bottom dotted lines. 
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FIG. 3. AE vs. for data in the —>■ D'^ D~ analysis. The signal region is indicated by a 
soUd box. The sideband region Hes above the top and below the bottom dotted lines. 
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